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CHAPTER 


INTRODUCTION 


e dynamic  beh 


reactivity  a 


a> 


reactor  to  measure  the  dynamic  parameters.  These  inves- 
tigations, while  giving  insight  into  the  fundamental 
nature  of  multiplying  systems  and  information  on  the 


Moots  a)  (8)  was  among  the  first  to  relate  the 


reactor  noise  spectrum  which 


of  the  effective  delayed  neutron  fraction,^}.  to  the 
prompt  neutron  lifetime,  £ . This  technique  has  become 


almost  standard  02)  • 


not  definable.  This  difficulty 


Ol)  develops 
iccording  to  th 


indie 


of  Rice  (13).  Rece 


»e.  Third,  the  autocorrelation  or  pile  noi 
method  does  not  yield  any  direct  informatl 


5xl0”4)  detector  efficiency  i 
physically  impossible  to  plac 


echniques.  Balcomb  (1£)  and  Rajagopal  (1£)  invest 


ivity  inputs . Both 
t impulse  response  f 


Is  giv 


will  reveal  something  about  the  nature  of  the  coupling 


regions.  The  cross-power  spectrum  between  the  outputs 
of  the  two  regions  Is  a complex  quantity  and  Its  phase 


disturbance  between  the  regions.  Its  amplitude  part 
as  well  as  the  phase  part  is  found  to  be  dependent  on 
the  value  of  the  multiplication  factor  of  the  undls- 


region  1,  and  (b)  no  external  Input  to  either  region. 


amplitude  mode 


CHAPTER 


THEORY 


and  their  Inputs  and  outputs 


as  to  be  directly  applicable  to  the  experiment.  This 


by 


jpled 


and  transfer  functions  together.  After  certain  simpli- 
fying assumptions  are  made,  the  equations  are  examined 


Cross-Power  Spectra 

Consider  a two  region  coupled  system  as  shown 

own  unit  Impulse  response  function  hx(t)  and  h2(t),  and 
outputs  given  by  fx (t)  and  f2 (t) . The  regions  are 
coupled  to  one  another  because  a certain  fraction,  A,, 
of  the  output  of  region  1 is  received  as  an  input  to 
region  2 after  a transit  time,  9 , and  vice  versa. 

Region  1 is  also  subjected  to  a random  stationary 

external  input,  i(t). 

of  region  1 and  region  2,  012(T),  is  defined  as  (19). 


0i2<r> 


lim  if' 

T-voo  J f].<t)f2(t+7)dt  . 


(1) 


en  specified 


where  T 


dlspla 


independe 


Figure  1.  Two  Region  System  Diagram 
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so  the  convolution  integral  between  the  input,  ou 
and  unit  impulse  response  function  is  given  as  (19) , 


flit) 


i 


’i<A> 


Input(t-A)dA  . 


(2) 


is  the  sum  of  the  external  input,  l(t),  and  the  frac- 
tion, A2,  of  the  output  of  region  2,  at  a previous  time 
corresponding  to  the  transit  time,  6 , between  the 


Input (t)  - i(t)  + A2f  2 (t  - $ ) . (3) 

Equation  (2)  is  now  written  as 

fi<t)  - hit A>  [i<t-A>  + *2t2<t-A-0)|  <jA, 

J-CO  1 J (4) 


A2f2(t-A-0)|  dAf2<t+r: 


(5) 


aquation  (5)  Is 


put  into 


*u<r» 


L(t-A>f2(t+T  )dtdA 

(6) 

f2(t-A-0!f2<t+T  IdtdA.  . 


the  external  input  and  the  output  of  region  2,  while  the 


-f_  hi<A></>12<r+A>'JA. 


integral  between  the  input,  output,  and  unit  impulse 


[(A)  Input (t-  A)d A 


(8) 


f2(t) 


-c 


In  region  2 the  only  input  is  that  fraction. 

*1.  of  the  output  of  region  1 which  reaches  region  2 
after  the  transit  time, 0 , between  the  regions.  Thus, 

input(t)  - Aifi(t-g).  (9) 

Equation  (8)  is  written  as 


f2(t) 


^(Aifxft-A-fltdA 


(10) 


which,  when  substituted  into  equation  (1),  gives, 


012<r>=. 


jjAJfiit+r-A-^idAdt. 


Once  again.  Inversion  of  the  order  of  integration  will 
yield  a more  convenient  grouping, 

fiitjfiit+r-A-gJdtdA . 

(12) 


in  which  the  inner  Integral  is  rec 
autocorrelation  function  of  the  output  of  region  1.  Equa- 


tion (12)  can  now 
with  equation  (7) 
output*  of  region 


f oo 

0 12  ( r ) ” *iy  h2(A)0u(T-A-0)dA.  (13) 

An  obvious  simplification  is  the  case  of  no 
external  input  to  the  two  region  system.  Bquation  (13)  is 


Investigators  (1£)  (17)  choose  to  measure  cross- 


Many 


15 


digital  computer.  If  the  cross-power  and  power  spectra 


eliminated.  The  equipment  available  to  the  author  alBO 


power  spectra,  aquations  (7)  and  (13)  must  then  be 


by  multiplying  both  sides  of  the  equation  by 

where  J = V -1  , U is  the  angular  frequency  (radians 


iirf_  feine‘,u'*r 

*2  r °°  r °° 

e~‘UT,'T] m »iai0ja(T.A*em. 


i£12<u)  - 

dial T' 


IL’e"1* 


isform  pair  defined  by  (20) 


w[.~  ^ 


iT)e~3UT  *t 


i2iU)ejUT dU  ■ 


(17) 


3ft  side  of  equation  (15)  is  then  by  defini- 


• 12 ((*/).  In  the  first  double  Integral  on 


x-  r ♦ A . 


(18) 


(1„ 


Equation  (19)  becomes  a product  of  two 


KAMA]  £ jffJ  012(x)g‘J^xdxj  . 

(20) 


the  right  of  equation  (15)  the  change  of  variable 


y * T * A + 9 . (23) 


which  will  give  a product  of  two 


as  in  the  first  integral 
single  integrals, 


complex  comji 
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of  the  transfer  function  in  region  1,  and  the  second  factor 
in  brackets  is  the  power  spectrum  of  the  region  2 output. 
Thus  equation  (15)  may  finally  be  written, 

<|>12<U>  - Hj(L/)^i2((J)  ♦ A2ejL^H*(U)<f22<L/>.  (25) 

Equation  (13)  may  be  manipulated  in  the  same  manner 
as  was  equation  (7).  Taking  the  Fourier  transform  of  both 
sides  of  equation  (13)  gives, 

<t>i2lT)e']uT  aT 

e"JLyrdT/<D  h2a,0ii(r-A-0>dA. 

(26) 

(26)  to. 


z = T - A - 9 (27) 

will  give,  after  a separation  of  variables. 


(28) 


(29) 

(29) 


$„iu)  -*le'J‘j9vs(U)(|unji  . 

fer  functions  in  the  set  of  equations  (25)  and  (29) 
in  the  set  of  equations  (30)  and  (31) . The  reasons 


(31) 


symbolizing 


10  external  input  and  the 


4>u«ji  - ■1(ui§„hji  ♦ c 

“ Al6  iU^  (; 

Two  Region  Transfer  Funetiona 
The  two  transfer  functions,  H VIU) . and  H 


subcrltlcal  reactor  with  a neutron  leakage 
from  the  other  slab.  Within  each  region, 
space  independent  and  a one  group  bare  rea. 
assumed.  Cohn  Ql)  has  commented  that  the 
bare  reactor  model  Is  still  used  for  pract: 
kinetics  work,  even  though  it  has  been  dis< 


-ally  all 
irded  as 


8 1/v  lifetime . 


The  kinetics  equations  for  region  1 may  then  be 
written  (18) , 


I A±ctl 


^ikini 


Aicii 


n density  in  region  1 


ft-  l fit 


s the  concentration  of  precursors  of 
-th  group  of  delayed  neutrons  in  region  1 


caused  by  the  leakage 
-s  term  may  be  repre- 
► neutron  density  in 


the  Interaction 
it  from  region  2. 
sented  by  some  fraction,  a2,  < 
region  2 at  some  previous  time  corresponding 


(36) 
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desirable  to  retain  at  least  the  first  few  terms  of  the 
expansion  of  equation  (37) . 

The  kineticB  equations  for  region  1 are  then, 


•H  [yi-lH 
at  2 


A2n2  *2  6 dn  *20  2 d2n; 
i l dt  + 2 £ dt2  ' 


(38) 


£ 


(39) 


A similar  set  of  equations  exists  for  region  2.  in  the 
experiment,  a random  reactivity  input  was  given  to  region 
1 only,  while  *2  remained  constant.  This  fluctuation  in 
the  multiplication  factor  in  region  1 will  give  rise  to 
fluctuations  in  the  neutron  density  and  delayed  neutron 
precursor  concentrations  in  both  regions.  There  may  also 
be  fluctuations  in  the  external  neutron  source.  It  is 

to  linearise  these  quantities  as  follows: 


+ <5»i 

♦ 

* <5  cii 


“ n20  + 
“ C120  + 
“ S10  + 


<5»: 


convenient  then 


equations 


When  these  definitions  are  substituted  into  the 
of  region  1,  (38)  and  (39). 


d_  , [(*10+ 6 *!>(!-,#)-!]  (n10+(5ni) 

dt  10  1 i 

+ I^i<ciio+(5cii>  +sio+(5si  (41) 

J(n20+6n2)-^(n20+(5n2)+-^-^j(n20+(5„2,+  . . .j 

(42) 

JL,„  .,S_  . AO'lO^KiHn^dn,)  . 

at'  no  wcu>  y Vcil0+(5cu) 

The  suns  of  the  steady  state  portions  of  these  equations 
are  equal  to  zero,  and  the  equations  are  linearized  by 
setting  the  product.  6k fin.  equal  to  zero  since  both 
quantities  will  be  small.  The  equations  for  region  1 are 


d<5cji  ^1n1Q(Sx1  A*10  <5»,  , <• 

dt  i * a - Ai°ci 


conditions, 


i the  Laplace  transform,  using  zero  initial 
is  taken  of  both  sides  of  both  equations, 


sAh  - a-l)niAlti  , [»i»-l>-i]  A»,  r’,  . 

i g g +4-iu 

+ A<^1  + jf[A»2  - 0-A»  2 + 4"  -*A»2  - •••] 

^1n1Aic  fl.k.All. 

8 Acil  “ 7) + n ' ^i^Cii,  (46) 


where  s is  the  Laplace  transform  variable,  the  capital 
letters  denote  the  Laplace  transform  of  the  (5  quanti- 
ties (1  - e . , A is  the  Laplace  transform  of  6 02),  and 
the  small  letters  are  now  understood  to  represent  steady 
state  values.  Equation  (46)  may  be  solved  for  A C±1- 
When  this  is  done  and  A substituted  into  equation 
(45)  and  terms  collected. 


If  the 


multiplic 


reactivity  ti 


while  in  the  latter,  fluctuation  in  the  multiplication 


apply.  The 


the  ratio  of  the  output  of  slab  2 
1 when  both  are  in  Laplace  no 


When  a random  reactivity  input  is  given  to  region 
e desired  transfer  function  of  that  region  is  given 

Auj/n, 

H.tl.n  . i i 


where  A N1/n1  is  the  normalized  output,  and  A is  the 
input.  This  ratio  is  obtained  from  equation  (47)  after 
A la  eliminated  by  substitution  from  equation  (49). 
With  the  assumption  of  no  external  neutron  source,  solv- 
ing equation  (49)  for  A N2  gives, 


An2  » 


, 02*‘ 


>a(l-^)-l  h,  y 

l *9,  i s + A, 


Substitution  of  this  expression  for  A N2  into  equation 
(47)  gives,  after  some  algebraic  manipulations  and 
setting  s - J U and  Q - x - j(jO  4-  U202  . 

H,(L/>  = 


(48). 


u-j3i- 


*r? 


Ml. 


while  Y 2^(jJ ) Is  formed  directly  from  equation  (49),  as 
was  B 2IU). 


V2(U) 


Hi  AAj 
£ i ju+Aj. 


Reduction  of  Equations 


only  the  higher  frequencies  ( U > 1)  are 


1 - kl<l  - ^3  ) 

B1  = l ' ,S9) 


11  -P''l 

. JU  - ^<tl  e-“ue 


B2iU)  - Y (U)  - 


VI 

Bx  -f  jU 


(30),  (31),  (32),  2 
input  is  given  to  i 


§atv>  - 1 


e 


$12(U)  - — 


4]e-^fu 


, JU  - K»a^Jle'a;IL'e  , 


<i>12<(*/) 


[Vtl  e $ii(wi 


■dependent  of  direction.  Further  amplification 
d if  the  slabs  are  assumed  identical  except 
multiplication  factors  which  may  be  different, 
ptlon  results  in  the  condition 


f, H>)e-JU»  $n(U) 


pression  for  $u(t/)  when  there  is  no  external 
vlty  input  to  region  1 is  identical  to  equation 
wer  spectrum,  is  a real  quantity  but  t 

power  spectrum  is  a complex  quantity  which  may  t 
terms  of  its  phase  and  amplitude  componen 
complex  quantity  Rl2((J)  is  defined  as 

$la  iU) 


12 IU)  « 


nplltude,/R12/. 


and  its  phase,  p(R12),  parts,  as 


(70) 


p(r12) 


-Ucob9u  - B2sin0L/-| 

B2cob0U  - Usln0(J  J 


(71) 


Equation  (70)  shows  that  the  amplitude  portion  of  R12  is 
Independent  of  9 and  that  the  shape  of  / R12  / depends 


region  not  disturbed  by  the  random  reactivity  input. 
The  phase  portion  of  R12(U),  however,  is  strongly  de- 

9U  becomes  large.  Figures  2 and  3 show  plots  of  the 


of  the 


malized  amplitude  plot  (Figure  2)  is  made  for  several 
values  of  B:  100,  200,  500,  and  1000.  These  values 
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saaj&ap  '(/*))' 


ot6uv  aseqa 


slabs.  The  phase  plot  (Figure  3) 


indicates  the  dependence  of  the  complex  quantity  R12(L/) 
on  this  transit  time,  0 ■ For  B-500,  the  phase  plots 
for  three  different  values  of  9 are  shown:  first,  for 

0 equal  to  zero;  second,  for  0 corresponding  to  the' 


velocity;  and  third,  for  9 correspond- 


travel  to  the  other  side  with  the  velocity  of  a thermal 


Appendix  A.  The  phase 


part  of  R12  (l_J)  is  also  given  in  Figure  3 for  B values 
of  100,  200,  and  1000,  with  9 corresponding  to  the  ther- 


the  quantity  >*12 ((*/)  provides  a rather  sensitive  test 
for  the  validity  of  assumptions  about  the  nature  of  9 . 


mlnation  of  the  parameter,  B,  which  yields  information 


ch  the  multiplica tl 
Ml,  Baldwin  (H)  t 


Thus  a determination  of  B yields  information  about  the 


apllng  of  the 


CHAPTER 


slf  well  to  digital 


Results 


of  auto  or  crosscorrelation  measurements  are  in  terms  c 
the  unit  impulse  response  function,  while  the  quantity 


0 12*7")  -J  $12<L/>  e Jtjr  dU 


The  definition  of  the  crosscorrelation  funct 


hi ' 

T 

be  substituted  into  (17) 

■ -/  ’*> 


^(tjfjit+rjdt,  a) 
may  be  substituted  into  (17)  to  give, 


:(U)elUTdU 


(73) 


f2(t) 

slab. 


course,  be  performed  over  infinite  time  or  over  negative 
times,  so  a finite  average  is  used,  resulting  in  the 


approximation. 


§12iu)eiL/T du  . 


"T { mji  ■ 

(75) 


’l  'l'°  Ju'i|“ri  L “ - $u“J1ieIU*r  Au, . 


-r(  [fi,«]u«-$nlUil4ul  . 


b the  mean  square  value  of 


b expression  is  solved  i 


r T 

I [■' 


tAu, 


Figure  4 indicates  the  sequence  of  operations  necessary 
to  solve  equation  (78)  experimentally. 

The  cross-power  spectrum  is  a complex  quantity 
and  may  be  expressed  in  terms  of  its  real  parts,  Re12, 
and  its  imaginary  parts,  lm12.  Use  is  also  made  of 
Euler's  relation. 

p J LJjT 

C = co  bUJ  + J sinter  , (79) 


fi(t)  ] f2<t+n  1 dt  = 

Jo  1 J u1 

(■  (80) 

[Re12cosUir-im12sinUir. tjlim^cos  A 


fift) 


$11iu1k4u1 


),  equation  (80)  reduces 


r^'  ]u‘ 


The  integral  on  the  left  side  of  this  equation  is  a 
real  quantity  so  that  when  the  real  parts  are  equated, 


tAl/x 


Returning  to  equation  (80) , one  sets  T equal  to 

2U1  ' 

Since  both  f^ (t)  and  fj (t)  have  been  put  through 
identical  band  pass  filters  and  are  both  essentially  at 
the  same  frequency,  ij  y 


Thus  equation  (80)  becomes, 


(83) 


[-Im12+JRe12  j A Uv 


(84) 


§ 

I 

I 

I 


UFTR 


52 


10«6' 


fllte 


tx(t) 


U1 


lag  in  a low  pass  filter  and 
given  a phase  lead  in  a high  pass 


filter  such  that  there  was  a 


it  and  the  lagging  fx(t)| 


shifting  networks  were  tore  complex  than  the  simple 
filters  indicated  in  Figure  9 because  each  filter  fed 


/Vt 


v 


(R  + K)2  + rVuV 


(86) 


_ bkUc  1 


K 


V (R  + K)2  + RVUV 


(89) 


where  K Is  the  value  of  the  input  resistor  to  the 
operational  amplifier.  Since  the  signals  were  passed 
through  the  band  pass  filters  and  were  at  essentially  a 


give  the  desired  90  degree  phase  shift  at  each 


a RC  combina 
high  pass  an 


-h  that  the  amplitude,  j hj,  of  the 


57 

multiplication,  however,  the  signals  were  passed 
high  pass  filters  to  remove  any  dc  components 
frequency  drifts  introduced  because  of  the  high 
the  system.  The  nominal  1.0  microfarad  capaci- 
1 10  megohm  resistors  used  were  actually  1.09 


The  separate  and  combined 


applied  to  the  real  and  imaginary  parts  of  the  cross- 
power spectrum  are  given  in  Appendix  B. 


Multiplication  was  performed  after  the  sig 
assed  through  the  high  pass  filters  and  the 


the  cross-power  spectrum.  If  the  gain  of  the  opera- 
tional amplifiers  between  the  band  pass  filters  and  the 
multipliers  was  M.  the  gain  of  the  integrators  (1/RC) 
was  P,  the  voltages  on  the  integrators  were  V for  the 


band  width  were  f1  and  Afx,  respectively,  then  the 
real  and  imaginary  parts  of  |>12<U)  were  given  by. 


The  measurement  of  the  quantity  R 12 (U/)  as 
defined  by  equation  (69)  requires  that  both  the 


cross-power  spectrum  between  both  outputs  and  the  power 
spectrum  of  the  output  of  the  externally  excited  region 
be  measured.  The  real  and  imaginary  parts  of  the  cross- 
power  spectrum  were  measured  by  use  of  equations  (91) 
and  (92).  The  power  spectrum  of  region  1 was  measured 
with  the  same  computer  setup  as  was  the  cross-power 

quantity,  there  was  only  a voltage.  V.  on  the  Integra- 


T Af,(M)2P 


rhuB  the  real  and  imaginary  parts  of  R (f  ) were 
given  by, 


Re  [ R12  ( f 1 > ] 

In[R12<fl>  ] » — 


phlte 


where  g is  the  average  number  of  sign  changes  per 


The  power  spectrum  of  this  signal  was  then  made  white 
out  to  the  desired  frequency  by  a corresponding  increase 
in  the  count  rate  seen  by  the  GM  tube.  The  signal  was 
given  the  proper  amplitude  and  dc  level  on  the  analog 


The  reactivity  device  was  sufficiently  long  to 


CHAPTER 


EXPERIMENTS  AND  RESULTS 
This  chapter  contains  a description  of  the 
experiments  performed  both  for  check-out  and  calibration 


taken,  data  reduction,  and  the  confidence  limits  assigned 
to  the  results.  A detailed  description  of  the  procedure 


this  input.  The  error  analysis  associated  with  this 
calibration  is  used  with  the  data  obtained  from  the 


with  and  without  an 


led  both 
snted. 


r spectrum  analyzer 


the 


low  pat 


filter. 


the  output,  is  equal  to  the  product  of  the 

power  epectrum  of  the  input,  <^U(L/) , and  the  trans- 
fer function  of  the  system,  a(U)  09).  The  real  and 
imaginary  parts  of  the  transfer  function  may  be 
measured  by. 


<f>ii<u) 


(97) 


“(•H 


§ii<U) 


2 low  Pass  filter,  the  transfer  function  of 


the  filter  fed  a 0.1  megohm 
operational  amplifier  on  the 
amplitude  and  phase  parts  of 
G(L/) . were  calculated  from. 


(91) 


RESULTS  OF  RUN  L 


f Normalized 

Amplitude 

cycles  per  o£  GtU) 
second  decibels 


PofSG  (U)le 
degrees 


Normalized  Power 
Spectrum  of 

decibels 
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saoiBap  '(f))g  jo  st6uv  sse<u 
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3pn3TT*HY  P8ZTT*“zoK 


Cf/N 


(99) 


where  T ia  the  integration  time  and  A f is  the  frequency 
difference  between  half-power  points  in  the  band  pass 


■M 


which  gives  the  standard  deviation,  O' f.  of  a function 
f (x,y)  of  the  two  quantities  x and  y which  have  standard 
deviations.  0"x  and  0"y , respectively. 

The  above  expression  gives  the  standard  deviation 
of  the  phase  part  (given  by  the  quotient  VVvR)  aB, 


»»  V? 


the  amplitude  part  at  the  lower  frequencies  is  rather 
large.  The  standard  deviation  of  the  phase  part,  how- 
ever, is  of  a magnitude  comparable  with  the  size  of  the 


recording,  and  demodulation.  The  magnitude  of  the 
amplitude  scatter  le  about  ±2  decibels,  while  little 
difference  is  noted  in  the  phase  part.  In  all  other 


Run  P demonstrates  that  the  techniques  and 


RESULTS  OP  RUN  P 
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if  it  represents  the  end  of  correlation.  The  power 
spectrum  of  the  region  1 output  indicates  that  corre- 
lation would  end  much  beyond  25  or  30  cycles  per  second. 
At  30  cycles  per  second  the  power  spectrum  is  down  40 
decibels  from  the  amplitude  at  the  lower  frequencies. 
From  the  calibration  runs  it  would  be  expected  that 
correlation  would  not  continue  when  the  input  (which  is 
the  output  of  the  south  fuel  region  in  this  case)  power 


spectrum  drops  off  40  decibels.  There  was  no  drop-off 
of  the  spectrum  at  the  lower  frequencies  as  in  Run  R. 


ity  input  device.  The  magnitude  of  the  input  signal  was 
also  greater  than  before  resulting  in  the  lower  fre- 


quency portion  of  the  spectrum  being  50  decibels  above 
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The  amplitude  data  points  clearly  indicate  a 
loss  of  correlation  beyond  frequencies  of  20  cycles  per 
second.  The  cross-power  spectrum  between  the  ion  chant- 
ers should  be  identical  to  the  output  power  spectrum 
of  the  region  if  there  is  sufficient  signal  above  the 

amplitude  part  of  this  cross-power  spectrum  is  nor- 
malized by  the  output  power  spectrum  of  the  region  as 
determined  by  one  of  the  chambers,  a constant  amplitude 


t amplitude  is  the 


by  the  phase  part  data  points.  There  would  be  no  Dh 


tion  existed.  The  zero  phase  angle  line  is  plotted  i 
figure  25.  Beyond  20  cycles  per  second,  the  data 
scatter  indicates  a loss  of  correlation.  The  output 
power  spectrum  of  the  south  region  as  determined  by 
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corresponds  to  a multiplication  factor  of  0.874.  criti- 
indicate  that  the  reactor  would  be  critical  with  a one 


ing  of  3225  grams  of  uranium-235,  or  1612. S grams  per  slab, 


y strongly  coupled,-  t 


worth  about  3S0  grams  of  uranium-235.  As  much  as  1950 


in  a one  slab  configuration  gave  not  only  an  estimate 


s but  also  an  estimate  of 
given  mass  of  fuel.  For 


of  hay  25,  1961,  0.87*0.08,  loading  of  March  21,  1962. 


on  factor  of 


by  the  cross- 
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APPENDIX 


CALCULATION  OP  TRANSIT  TIME  BETWEEN  SLABS 


that  tha  neutron  wave  originates  in  the  center  of  one 
fuel  slab  and  travels  to  the  center  of  the  opposite 
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Thompson  (29)  gives  the  following  values  at 
thermal  energies  for  the  fuel  and  graphite  regions  of 


Graphite 


These  results  are  given  for  a fuel  slab  consisting  of  a 
homogeneous  mixture  of  water,  aluminum,  and  uranium. 

It  is  also  assumed  that  the  reactor  is  loaded  with  3500 
grams  of  the  20  per  cent  enriched  fuel. 

For  the  reactor  loaded  with  3100  grams  of  the 
fully  enriched  fuel.  £ for  the  fuel  section  is  0.0723 
centimeters  . The  diffusion  coefficient  and  diffusion 
length  of  the  fuel  region  will  also  change  from  the 
values  given  for  the  20  per  cent  enriched  fuel.  These 


changes  are  small  enough,  however,  and  the  lag  time  ] 
the  fuel  is  a small  enough  fraction  of  the  total  lag 
time,  that  the  values  for  the  20  per  cent  enriched  ft 
may  be  used  for  the  fully  enriched  fuel  also. 

The  dimensions  a and  b are  taken  as  264  and  1 
centimeters.  These  dimensions  correspond  to  the  com- 
plete diffusion  interface  and  include 
ite  surrounding  the  fuel  slab.  The  n« 
taken  as  2200  meters  per  second. 


>f  the  graph- 
velocity  is 
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i the  fuel  section  then, 


6’  - ^ i . 


Since  the  wave  must  travel  through  halves  of  two 
fuel  sections,  the  lag  time  through  one  full  fuel  sec- 
tion is  calculated.  This  lag  time  may  be  found  by 
dividing  the  thickness  of  the  fuel  region,  13.34  centi- 
meters, by  the  wave  velocity  as  determined  by  equation 
(104).  Substitution  of  the  fuel  region  value  gives. 


L.78xlo“)%U2]  - 1.78x10“  i 


is  useful  in  the  evaluation  of  0£u<jl  for  small  values 
oiU.  me  value  of  the  lag  time  in  the  fuel  at  various 
frequencies  is  given  in  Table  8. 


Ill 


TABLE  8 

VALUE  OF  TRANSIT  TIME  BETWEEN  SLABS 


cycles/ second 


Lag  Time  In  Milliseconds 
Fuel  Graphite  Total 


10 

20 

50 

100 

200 


2.60 

2.56 

2.20 


The  value  of  (5  2 in  the  graphite 


142 J I 51.8J 

(52  - 0.001012  . 

The  thickness  of  the  graphite  section  separating  the  two 
fuel  regions  is  30.48  centimeters.  When  this  distance 
is  divided  by  the  velocity  of  a neutron  wave  in  graphite 
as  determined  by  equation  (104) . the  quotient  is  the 
lag  time  in  the  graphite. 


[(2.04x102)2+L/2]  -2.04xl02j 


the  lag  time  in  the  graphite 
e also  given  in  Table  8. 
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The  R and  C values  at  each  frequency  set 

Polystyrene  capacitors  were  used  in  values  of  1, 
0.1,  and  0.01  microfarads.  The  R values  were  set  on  10 
turn  0.1  megohm  potentiometers.  Any  additional  resist- 
ance was  provided  by  adding  precision  resistors  in 
series  with  the  potentiometers.  These  tl%  resistors 
were  available  in  values  of  0.1,  0.2,  0.5,  and  1.0 
megohms.  The  potentiometers  were  used  as  variable  re- 
sistors and  their  set  points  are  also  given  in  Table  9. 

The  attenuation  of  the  signals  through  the  phase 
changing  filcers  was  calculated  from  equation  (110)  for 
the  two  values  of  input  resistors  used,  0.1  and  1.0 
megohms,  when  this  attenuation  was  normalized  to  one 
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£or  the  complete  isolation  of  the  filter,  i.e.  for  K 
very  large,  the  reciprocal  squared  of  the  attenuation 
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